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INTRODUCTION 

Mineralogical  applications  of  reaction  calorimetry  may  involve  either 
(al  direct  combination  and  transformation  studies,  or  (b)  solution  calorimetry. 
Among  these  approaches  the  latter  is  the  more  generally  useful.  While  solution 
calorimetry  using  acid  solvents  at  or  near  room  temperature  has  proved  useful 
in  some  cases,  it  has  very  serious  limitations  when  the  minerals  involved  are 
dense  and  refractory  This  is  the  principal  domain  of  high-temperature  solution 
calorimetry  and  of  the  modern  heat  flux  calorimeter. 

One  of  the  crucial  problems  in  high-temperature  solution  calorimetry  is 
the  choice  of  a suitable  solvent,  and  some  criteria  used  in  the  selection  of  oxide 
melt  solvents  for  work  on  refractory  minerals  are  outlined.  Examples  of  applica- 
tions of  reaction  calorimetry  to  mineral  systems  include  studies  of  (a)  enthalpies 
of  formation  from  component  oxides,  (b|  phase  transformations,  (c)  disorder 
problems,  and  (d)  mineralogical  reactions  and  equilibria. 

For  reacting  system  at  fixed  composition  or  for  a simple,  first  order  phase 
transformation,  we  have  the  well-known  thermodynamic  relation 


AGT  = AC"  + ^ MrdP 

rr 

= AH?-TAS"r+\ 


A VTdP 


At  equilibrium  ACr  = 0.  Thus  we  have  also 


i: 


AlfdP  = - AGj 


(I) 


(21 


In  these  expressions,  AG"  is  the  standard  Gibbs  energy  change  at  temperature 
7.  A H?  is  the  standard  enthalpy  and  AS"  is  the  standard  entropy  change; 
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A l , is  the  volume  change  associated  with  the  reaction.  Due  to  compressibility 
Al  t will,  in  general,  be  different  from  A l7?  the  volume  change  at  one  atmosphere 
pressure.  However,  as  long  as  we  are  dealing  with  condensed  phases  and 
moderate  pressures,  we  have 

Alr  % AV? 


and 


, ag? 
AF° 


(31 


Thus,  if  AG  ” and  Al , are  known.  may  be  calculated  or  estimated. 
The  Clapeyron  equation 


similarly  provides  a link  between  the  temperature  derivative  of  the  equilibrium 
pressure  and  the  change  in  entropy. 

The  implication  of  these  relationships  is  clear.  If  we  have  information  on 
AH0.  AS",  uml  Al  " for  a given  reaction  or  phase  transformation,  i.e.  if  we 
have  thermodynamic  data  valid  at  one  atmosphere  pressure,  we  can  readily 
estimate  the  pressure  temperature  equilibrium  curve.  If  we  also  know  Al 
as  a function  of  temperature  and  pressure,  we  can  obtain  more  precise  values. 

Our  principal  interest  in  the  present  paper  will  be  focused  on  the  purely 
thermal  quantities  AS0  and  A H°.  For  completely  ordered  phases  S°  and 
AS0  in  principle  may  be  derived  from  low  temperature  heat  capacity  data, 
via  the  third  law 

Where  the  relevant  heat  capacity  data  are  lacking,  approximate  values  of 
S°  for  condensed  phases  at  high  temperatures  often  can  be  estimated  However, 
complications  arise  where  the  substances  involved  are  disordered.  Under  these 
circumstances  configurational  entropies  sometimes  can  be  estimated  from 
detailed  structural  information. 

It  is  more  difficult  to  make  realistic  predictions  or  estimates  of  A Hn.  which 
usually  must  be  obtained  from  data  based  on  reaction  calorimetry. 


SOLUTION  CALORIMETRY  OF  MINERALS 

In  the  years  immediately  following  World  War  11,  considerable  progress 
was  made  in  the  general  area  of  mineral  thermochemistry.  This  progress  was 
sparked  by  the  development  by  Torgeson  and  Sahama  ( 1 948)  of  precision 
hydrofluoric  acid  solution  calorimetry.  The  thermodynamic  basis  of  this 
method  is  very  straightforward  . 
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A + B + Solvent  = Solution : AW , (7a) 

AB  + Solvent  = Solution:  AW^  (7b) 

A + B = AB;  AW"  = AW,  — AW2  (7c) 

Thus  the  standard  enthalpy  change  for  a given  reaction  (here  A + B = AB) 
may  be  obtained  from  the  difference  between  two  measured  enthalpies  of 
solution.  Following  the  work  of  Torgeson  and  Sahama.  hydrofluoric  acid 
solution  calorimetry  was  adopted  by  a number  of  investigators  in  the  silicate 
field,  notably  in  studies  of  feldspars  by  Kelley  el  ai.  ( 1953).  Kracek  and  Neuvonen 
(1952).  and  others.  More  recent  work  based  on  this  method  has  been  carried 
out.  e g.  by  Waldbaum  (1966). 

It  is  believed  that  the  first  purely  thermochemical  calculation  of  the  P T 
curve  for  an  important  silicate  equilibrium  was  that  of  Kracek.  Neuvonen 
and  Burley  (1951)  on  the  stability  of  jadeite.  This  work  preceded  by  several 
years  the  experimental  equilibrium  work  of  Robertson,  Birch  and  MacDonald 
(1957).  and  of  Birch  and  Le  Conte  (I960).  While  the  agreement  between  the 
calculated  and  the  experimental  values  was  far  from  perfect,  these  studies 
served  to  illustrate  both  the  potential  of  the  thermochemical  approach,  and 
also  some  of  its  obvious  limitations.  Thus.  Kracek.  Neuvonen  and  Burley 
found  that  jadeite.  as  most  other  dense  minerals,  dissolves  only  with 
considerable  difficulty  in  hydrofluoric  acid,  even  at  75  C.  Their  results  also 
indicated  that  the  observed  enthalpies  of  solution  to  some  extent  reflected  the 
past  treatment  of  the  jadeite  samples,  i.e.  whether  the  jadeite  powders  had 
been  ground  in  an  agate  or  a mullite  mortar.  Another  problem,  w hich  is  common 
to  all  work  on  hydrofluoric  acid  solution  calorimetry,  is  the  fact  that  the  observ- 
ed enthalpies  of  solution  are  large  numerical  quantities  compared  to  the 
enthalpy  changes  in  typical  mineral  reactions  and  phase  transformations.  As 
a result  of  these  various  problems  it  soon  became  recognized  that  HF  solution 
calorimetry,  although  clearly  useful  for  certain  types  of  mineral  equilibrium 
studies,  has  a somewhat  limited  range  of  application  A new  and  more  powerful 
approach  was  provided  by  oxide  melt  solution  calorimetry. 

In  this  method,  which  originated  with  the  work  of  Yokokawa  and  Kleppa 
(1964.  1965).  the  aqueous  a>id  solvent  of  conventional,  room  temperature 
solution  calorimetry,  is  replaced  by  a high-temperature  oxide  melt  solvent. 
The  technique  has  certain  obvious  advantages  compared  to  the  room  tempera- 
ture methods: 

( 1 ) It  allows  thermochemical  study  of  a wider  range  of  oxidic  substances,  and. 
particularly  of  dense  and  refractory  minerals. 

(2)  The  observed  enthalpies  of  solution  are  numerically  much  smaller  than 
in  acid  solution  work.  Without  loss  in  the  precision  of  the  enthalpies  of  reaction 
this  permits  the  use  of  heat  flux  calorimetry. 

(3)  It  is  possible  to  work  with  much  smaller  samples  than  in  conventional 
acid  solution  calorimetry.  With  samples  of  the  order  of  25  mg,  it  is  feasible  to 
combine  solution  calorimetry  with  high  pressure-high  temperature  synthesis 
and  equilibrium  work. 
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However,  there  are  also  disadvantages.  For  example,  it  is  not  convenient  to 
use  oxide  melt  solution  calorimetry  for  the  study  of  the  hydrous  minerals 
which  frequently  participate  in  low-grade  metamorphic  reactions. 

REACTION  CALORIMETRY  AT  HIGH  TEMPERATURE 

Two  basically  different  ideas  have  provided  the  basis  for  design  of  calorimeters 
operated  at  high  temperature.  The  idea  behind  the  adiabatic  and  ‘comentionar 
lor  quasi-adiabatic)  designs  has  been  to  reduce  as  far  as  possible  (ideally  to 
eliminate)  the  heat  (lux  (‘heat  leak  ) between  the  calorimeter  and  the 
surroundings. 

There  are  several  different  circumstances  which  tend  to  limit  the  effectiveness 
both  of  adiabatic  and  quasi-adiabatic  calorimeters  at  high  temperature 
The  most  serious  undoubtedly  are  related  to  the  increasing  importance  of 
radiation  in  heat  transfer  as  the  temperature  is  increased. 

In  the  design  of  heat  flux  calorimeters  no  attempt  is  made  to  reduce  the  heat 
leak.  Instead  one  tries  to  fix  the  design  parameters  ir.  such  a way  that  the  inte- 
grated heat  flux  can  be  measured  reproducibly  and  conveniently.  Figure  I 
illustrates  typical  time  temperature  curves  for  adiabatic,  conventional  and 
heat  flux  calorimeters.  In  our  own  work  we  have  found  it  convenient  to  design 
our  heat  flux  calorimeters  in  such  a w ay  that  the  time  constant  of  the  calorimeter 
is  of  the  order  of  3 minutes.  Under  these  conditions  a temperature  change 
which  occurs  in  the  calorimeter  at  time  0,  will  have  decayed  to  less  than  1 1000 


I igurc  I Schematic  time  temperature  curves  for 
adiahatic.  conventional  (i.c.  ‘quasi-adiabatic')  and 
micro  li  e.  'heal  flux')  calorimeters 


I 
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of  Us  maximum  value  in  about  one-half  hour.  This  usually  is  a convenient 
lime  scale  for  solution  calorimetry. 

Under  favourable  circumstances  Ifor  each)  the  heat  flux  calorimeter  cannot 
compete  with  the  adiabatic  calorimeter  with  respect  to  experimental  precision. 
In  fact,  we  have  found  it  very  difficult  to  achieve  a precision  better  than  about  + 
0.5  per  cent  in  heat  flux  calorimetry  carried  out  at  any  temperature.  However, 
while  the  difficulties  associated  with  adiabatic  and  conventional  calorimetry 
increase  dramatically  with  temperature,  those  of  heat  flux  calorimetry  increase 
much  more  slowly  Also,  what  heat  flux  calorimetry  may  lose  in  precision, 
it  may  regain  in  simplicity,  convenience  and  reliability. 

In  solution  (as  distinct  from  bomb)  calorimetry,  the  calorimetric  process 
usually  is  not  instantaneous,  but  may  require  a significant  period  of  lime.  In 
order  to  monitor  the  progress  of  the  dissolution  of  the  solid  sample,  and  the 


figure  2 Schematic  diagram  of  Calvet-type  twin  microcalori- 
meter  for  tempcraiure  up  to  about  800  C.  A:  Main  heater; 
B lop  heater ;(  Heavy  shield.  D Nickel  block . £:('alorimeter. 
F Protection  lube;  Ci  Radiation  shields.  H:  Manipulation 
lube;  + - : Thermopile 
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completion  of  the  dissolution  process,  it  is  desirable  to  use  a heat  flux  calorimeter 
with  very  good  baseline  stability.  Such  an  apparatus  is.  e.g.  the  Calvet-type 
twin  micro-calorimeter  shown  in  Figure  2 (Kleppa.  19MI)  The  twin  construc- 
tion. with  two  essentially  identical  calorimeters  and  thermopiles  connected 
in  series  but  backed  against  each  other,  is  particularly  useful  in  reducing  the 
need  for  very  close  temperature  control.  This  can  be  a very  difficult  problem 
at  high  temperatures.  The  apparatus  illustrated  in  Figure  2 has  a calorimeter 
block  ('thermostat')  constructed  from  pure  nickel,  which  is  a suitable  material 
for  operation  in  air  at  temperatures  up  to  7- 800  C.  By  choosing  other  construc- 
tion materials  such  as  inconel,  stainless  steel,  nichrome.  Kanthal.  or  aluminium 
oxide,  the  upper  temperature  limit  of  the  apparatus  may  be  increased  to  that 
indicated  for  the  chosen  material.  As  the  temperature  of  operation  is  raised, 
problems  frequently  arise  in  connection  with  electrical  noise  in  the  thermopile 
system.  Careful  shielding  of  the  thermopile  system  is  indicated. 

Hollow 
quartz  tube 

Hollow  quortz 
manipulation  tube 

Quartz  shield 

Top  level  of 
colorimeter 

Quartz 
outer  liner 

Quartz  crucible 

Plotinum  sleeve 


Gold  or  Plotinum 
crucible 

Plotinum 
sample  cup 


2Pb0-B?0j 
oxide  melt 


Koowool  plug 


figure  3 Schematic  diagram  of  experi- 
mental arrangement  used  in  oxide  melt 
solution  calorimetry. 
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Various  experimental  arrangements  have  been  tried  within  the  calorimeter 
proper.  Our  own  preferred  set-up  is  shown  in  Figure  3.  Prior  to  the  solution 
experiment  the  powdered  sample  (150  250  mesh;  - 25  mg)  is  held  in  a shallow 
platinum  cup  above  the  surface  of  the  solvent.  Dissolution  of  the  sample  is 
initiated  by  dipping  the  cup  into  the  melt.  Stirring  is  manual  by  moving  the 
platinum  cup  up  and  down  several  times. 


Oxide  Melt  Solvents 

In  their  initial  work.  Yokokawa  and  Kleppa  (1964.  1965)  explored  the  poten- 
tial of  liquid  vanadium  pentoxide  at  a calorimetric  solvent  at  700  C.  The 
results  of  this  early  exploration  are  summarized  in  Figure  4.  which  gives  the 


f igure  4 Molar  integral  enthalpies  of  mixing  of 
liquid  V2()s  + solid  oxides  at  about  700  C (Reprinted 
with  permission  from  T Yokkoawa  and  <)  J Kleppa. 
I nor  a ( hem  . 4,  1X06 1 1%5).  Copyright  hi the  American 
( hemical  Society). 
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enthalpy  changes  associated  with  the  solution  of  some  simple  oxides  MO 
in  V,0<.  Note  in  particular  that  for  the  oxides  of  the  alkaline  earth  metals 
the  enthalpies  of  solution  change  systematically  in  the  sequence 

BeO  > MgO  > CaO  > SrO  > BaO 

with  small  positive  values  for  BeO  and  quite  large  negative  values  for  BaO. 
This  is  of  course  the  order  of  increasing  ionic  radius  of  the  cation  and  of  increas- 
ing basicity  of  the  solute  oxide,  w hich  is  consistent  w ith  the  fact  that  the  solvent 
V ,0,  is  a very  strongly  acidic  oxide.  These  early  results  immediately  suggest 
that  an  understanding  of  the  elements  of  acid  base  chemistry  of  melts  is 
essential  in  any  meaningful  discussion  of  possible  solvents  for  oxide  melt  solu- 
tion calorimetry. 

It  is  convenient  to  discuss  the  acid  base  character  of  oxides  and  oxidic 
melts  in  terms  of  the  Lux  Flood  approach,  which  is  analogous  to  the  well- 
known  Bronsted  Lowry  scheme  for  protonic  acids.  According  to  Lux  (1939) 
and  Flood  and  F0rland  (1947),  a suitable  measure  of  the  basicity  of  the  melt 
is  the  oxygen  ion  activity,  which  is  determined  by  (he  equilibrium 

Acid  + O 2 = Base. 

The  stronger  the  base,  the  higher  the  oxygen  ion  activity. 

For  the  vanadium  (V)-oxygen  system  we  may  in  this  scheme  write  a series  of 
possible  (but  not  necessarily  realized)  consecutive  acid  base  steps,  which 
involve  species  such  as  V**,  VO°.  V02,  V2Ov  VO,.  V207  *.  and  V04  3. 
Similarly  we  may  for  each  oxide  MO  consider  the  species  M * 2.  MO.  MOj 2, 
etc.  When  a basic  solute  MO  is  added  to  liquid  V2Ov  this  will  give  rise  to  the 
formation  of  vanadate  species  such  as  VO,,  V207  4,  and  V04  3 (or  polymeric 
forms  of  these).  To  a first  approximation  the  magnitude  and  sign  of  the  enthalpy 
change  associated  with  the  solution  process 

MO(.v)  + V,0,(/)  = M f 2 + Vanadate  (in  V20,) 

may  be  considered  a measure  of  the  tendency  of  this  acid  base  reaction  to 
proceed  from  left  to  right,  i.e.  of  the  relative  strength  of  V20,  and  M * 2 as 
acids  or  of  MO  and  vanadate  as  bases.  The  stronger  the  base  MO  (or  the  weaker 
the  corresponding  acid  M ' 2| the  more  negative  the  partial  enthalpy  of  solution 

(AH  M(>l- 

If  the  solute  oxide  MnO„  rs  very  acidic,  e g.  has  an  acidity  which  approaches 
that  of  (he  solvent  V2Ov  the  acid  base  character  of  the  dissolution  process 
in  large  measure  will  be  suppressed.  Under  these  conditions  the  enthalpy  of 
solution  no  longer  will  serve  as  a meaningful  measure  of  the  relative  acid  base 
strengths  of  the  solutes.  For  this  purpose  one  should  of  course  instead  use  an 
oxidic  solvent  of  basic  character,  such  as.  for  example,  liquid  PbO. 

The  vanadate  V20,  mixtures  which  result  from  the  dissolution  of  basic- 
oxides  in  V2(),  are  examples  of  buffered  acid  base  melts.  Since  most  mixed 
oxide  compounds  contain  two  or  more  components  of  different  acidity,  it 
usually  will  be  most  convenient  to  use  buffered  melts  as  calorimetric  solvent. 


This  also  has  the  practical  advantage  that  in  the  buffered  regions  the  enthalpies 
of  solution  of  acidic  or  basic  oxides  tend  to  vary  less  with  composition. 

We  shall  illustrate  these  problems  further  by  considering  the  binary  liquid 
system  PbO  B,0,.  At  700  8(X)  C,  which  is  a convenient  temperature  for 
oxide  melt  solution  calorimetry,  the  liquid  range  in  this  system  extends  from 
about  15  to  100  per  cent  B,0,  (Geller  and  Bunting,  1937).  However,  there  is  a 
miscibility  gap  in  the  high  B20,  region,  which  is  also  characterized  by  a very 
high  viscosity.  The  system  is  glassforming,  which  has  proved  an  advantage 
in  checking  that  oxide  samples  do  in  fact  go  into  solution. 

Holm  and  Kleppa  ( 1967)  published  the  results  of  a calorimetric  investigation 
of  this  system  carried  out  at  800  C.  In  this  work  the  partial  enthalpies  of 


PbO  B2Oj  melts  at  X(X>  C □ Rodigina.  Gomel'skii, 
and  Lugimna  (1961); By  Gibbs  Duhcm  inte- 
gration from  (Reprinted  with  permission  from 

J.  L.  Holm  and  O.  J Kleppa,  Inory  ('hem.. 6, 645(19671. 
Copyright  by  (he  American  Chemical  Society). 
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solution  {actually  AH  An)  of  PbOt.sland  B20,(/)  were  measured  at  compositions 
between  about  17  mole  per  cent  and  about  67  mole  per  cent  B203.  At  higher 
BjO,  contents  the  melts  were  too  viscous  to  allow  precise  calorimetry  by  the 
direct  mixing  approach. 

The  results  reported  by  Holm  and  Kleppa  are  shown  in  Figures  5 and  6. 
Since  the  partial  enthalpies  of  both  components  were  measured,  it  was  possible 
to  check  the  internal  consistency  of  the  result  by  means  of  the  Gibbs  Duhem 
equation.  The  results  of  this  check  are  indicated  as  broken  lines  in  the  figures. 
On  the  whole  there  is  reasonable  agreement  between  the  calculated  and  measur- 
ed data. 

The  data  suggest  interesting  changes  in  the  acidity  of  the  melt  as  the  composi- 


in  PbO  B.O,  melts  at  800  C By  Gibbs 
Duhem  integration  from  /.Pb<).  {Reprinted  with 
permission  from  J.  L.  Holm  and  O.  J.  Kleppa. 
Inory.  Chem.,  6,  645  (1967).  Copyright  by  the 
American  Chemical  Society). 
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lion  varies  from  basic  in  ihe  high  PbO  region  lo  acidic  in  the  high  B20,  range. 
We  interpreted  lhe  results  to  indicate  that  the  liquid  lead  borate  system  contains 
at  least  two  complex  borate  anions,  which  have  oxygen  boron  ratios  of  about 
3.5  4.0  and  1.X4  1.91.  respectively.  These  ratios  were  obtained  by  attempting 
to  locate  the  two  steeply  rising  parts  of  the  partial  enthalpy  versus  composition 
curves.  Between  these  two  compositions  the  curves  are  flatter,  which  presumably 
reflects  the  presence  of  both  of  the  two  borate  anions  in  comparable  amounts, 
i.e.  a buffer  region.  Most  of  our  early  calorimetric  work  was  based  on  the  use 
of  melts  with  a M0:B20,  ratio  of  3: 1.  However,  as  a result  of  the  work  by 
Holm  and  Kleppa  we  now  prefer  a somewhat  more  acid  melt  with  a ratio  of 


Ylincralogical  Applications 

The  mincralogical  applications  of  high  temperature  reaction  calorimetry 
fall  into  the  following  main  but  related  categories:  (a)  enthalpies  of  formation 
of  refractory  mixed  oxides,  (b)  phase  transformations,  (c)  disorder  problems, 
and  (d)  mineralogical  reactions  and  equilibria. 

(a)  Enthalpies  of  Formation 

The  thermodynamics  of  formation  of  refractory  mixed  oxides  from  its 
component  oxides  is  not  easily  studied  experimentally  since  this  area  is  almost 
inaccessible  to  conventional  room  temperature  calorimetric  methods.  During 


figure  7.  Enthalpies  formation  of  2 3 spinels  from  component  oxides  (Reprinted 
with  permission  from  A Navrotsky  and  O.  J Kleppa.  J Inorg.  Nucl.  t hem.  30, 
47lJ|l9f>8)  Pergamon  Press  i.td.) 
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the  past  10  15  years  some  progress  has  been  achieved  through  solid  gas 
equilibrium  methods,  and  in  particular  by  means  of  the  solid  oxide  electro- 
chemical cells  introduced  by  Kiukkola  and  Wagner  (1957).  However,  as  is  so 
commonly  the  case,  while  work  of  this  type  may  provide  reliable  data  on  the 
Gibbs  energy  of  formation,  it  generally  does  not  allow  an  unambiguous  sepa- 
ration of  AO,  into  its  appropriate  enthalpy  and  entropy  components.  Oxide 
melt  solution  calorimetry  has  represented  a new  ‘handle-  on  this  problem  A 
particularly  significant  example  is  the  investigation  by  Navrotsky  and  Kleppa 
1 1968)  of  the  thermodynamics  of  spinel  formation,  from  which  we  take  f igure  7, 
which  gives  enthalpy  of  formation  data  for  a senes  of  refractory  2 3 spinels 


(MAI204),  tungstates  IMWO,).  carbonates 
IMCOj).  sulphates  (MS04)  from  component 
oxides  (Reprinted  with  permission  from  A 
Navrotsky  and  O.  J.  Kleppa.  Inorq.  Chem.,  8. 
756(1969).  Copyright  by  the  American  Chemical 


Society). 
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taken  from  the  families  MAI2()4.  MCia;04.  MMn204.  and  MFe,04.  More 
recently  this  work  was  extended  by  Muller  and  Kleppa  (1973)  to  the  corres- 
ponding compounds  MC'r,04.  which  were  studied  by  solution  calorimetry 
at  9()0  C.  The  systematic  variation  of  A//f  with  the  cation  is  of  considerable 
interest.  However,  Figure  8 which  gives  the  enthalpies  of  formation  of  the 
compounds  MAI,04  compared  with  corresponding  data  for  MW()4.  MCO„ 
and  MS04  shows  that  this  variation  with  cation  is  nothing  specific.  The  gene- 
rally parallel  character  of  the  curves  is  apparent.  The  shift  in  A//,  to  more 
negative  values  in  the  series 

mai,o4  < MWC)4  < MCO,  < mso4, 

simply  reflects  the  increasing  acidity  of  the  second  oxide  in  the  sequence 
AFO,  < WO,  < C02  < SC)3. 


figure 9 P T stability  fields  for  the  coesile-stishovile  equili- 
brium and  for  the  metastable  equilibrium  quartz-stishovite 
calculated  from  thermodynamic  data  Sclar  ei  al.  refers  to  Sclar. 
Young.  Garrison,  and  Schwartz  (1962);  Wentorf  refers  to  Wenlorf 
(19621.  Stishov  refers  to  Stishov  (1963).  (Reprinted  with  permis- 
sion from  J L.  Holm.  O.  J Kleppa  and  E.  F.  Westrum.  Geochim. 

C osmot  him.  Ada,  31.  2289  (1967)  Pergamon  Press  Ltd.) 
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( h)  Phase  Transformations 


It  is  frequently  found  that  phase  transformations  in  solid  mineral  systems 
are  very  sluggish,  and  that  therefore  equilibrium  can  be  achieved  in  the  labora- 
tory only  at  quite  high  temperatures,  for  this  reason  quantitative  equilibrium 
studies  carried  out  at  temperature  may  be  very  difficult  and  at  times  even 
impossible.  However,  there  are  no  such  restrictions  on  the  thermodynamic 
calculations,  which  in  principle  will  apply  at  any  temperature. 

Examples  of  important  calculated  phase  equilibria  in  mineral  systems  are 
reported  by  Holm.  Kleppa  and  Westrum  (1967)  in  their  work  on  the  high 
pressure  phases  of  silica.  Figure  9 gives  the  thermodynamically  calculated 
P T curve  for  coesite  stishovite.  and  for  the  quartz  stishovite  transformation, 
along  with  some  static  equilibrium  data  for  the  former  equilibrium. 


0 200  400  600  800  1000 

Temp.'C 


figure  10  Comparison  of  calculated  P T curves  for  the  kyanite 
sillimanile  equilibrium  with  some  experimental  curves. 

Anderson  and  Kleppa  (1969): Holm  and  Kleppa 

(1966); Weill  ( 1966);  Richardson.  Bell,  and 

Gilbert  1 1968); Allhaus  ( 1967);  <[>  Newton  (1966)  (Repro- 

duced with  permission  from  P.  A.  M Anderson  and  O.  J Kleppa. 

Amer.  J.  Sci.,  267. 285  (1969).) 
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The  thermochemical  results  for  the  quartz  stishovite  transformation  are  of 
particular  interest,  since  this  equilibrium  occurs  in  the  stability  field  of  coesite 
and  cannot  be  studied  by  the  usual  static  equilibrium  methods.  However,  this 
transformation  takes  place  in  shock  wave  experiments  (e.g.  see.  Davis.  1972). 

Another  phase  transformation,  which  has  been  rather  extensively  studied 
by  high  temperature-high  pressure  methods,  and  for  which  good  thermochemi- 
cal  data  are  available  both  on  A//0  and  on  AS0,  is  the  kyanite  sillimanite 
equilibrium.  Figure  10  gives  a summary  of  some  of  the  P-T  data  for  this 
transformation  taken  from  Anderson  and  Kleppa  (1969).  Note  that  while 
there  is  now  reasonable  agreement  between  calculated  and  experimental 
equilibrium  pressures,  this  agreement  does  not  extend  to  the  dP  dT  slopes  which 
are  somewhat  higher  in  the  equilibrium  experiments  than  in  the  calculations. 
There  is  now  good  reason  to  believe  that  this  discrepancy  is  due  to  some  Al  Si 
disorder  in  the  sillimanite  structure,  which  makes  the  entropy  of  sillimanite 
somewhat  larger  than  given  by  the  third  law  value  (e.g.  see.  Holdaway.  1971). 


(c)  Disorder  Problem 

It  has  been  recognized  for  some  time  that  structural  disorder  is  an  important 
phenomenon  in  a number  of  mixed  phases  of  refractory  character  such  as. 
I e.g.  in  the  spinels,  mullite,  pyroxenes,  feldspars  As  an  example,  let  us  first 

' consider  the  disordering  problem  in  spinel.  Mg[Al]j04.  At  low  temperature 

the  equilibrium  state  of  this  mineral  has  all  the  magnesium  ions  located  on  the 
tetrahedral  interstitial  sites  in  the  essentially  close  packed  array  of  oxygen 
ions;  similarly,  the  aluminium  ions  are  located  on  octahedral  sites.  However, 
when  this  mineral  is  heated  above  800  C,  a disordering  transformation  occurs 
which  moves  a fraction  (x)  of  the  magnesium  ion  to  the  octahedral  sites  and  a 
corresponding  number  of  aluminum  ions  to  the  tetrahedral  sites  (Hafner  and 
Laves.  1961).  The  state  of  this  partially  disordered  spinel  may  be  described  as 

Mg,  ,AI,[AI}  tMg,]04. 

F igure  1 1 gives  a plot  of  the  configurational  entropy  associated  w ith  disorder 
in  a spinel,  as  a function  of  the  disorder  parameter  x.  Note  in  particular  the 
very  sharp  (infinite  slope)  rise  of  the  configurational  entropy  for  small  values 
of  x.  It  is  this  sharp  initial  rise  of  the  entropy  with  x which  makes  structural 
disorder  a very  general  problem  in  refractory  mixed  oxides  at  high  temperature. 

To  a first  approximation  it  may  be  assumed  that  the  endialpy  change  associat- 
ed with  the  disordering  process  in  spinels  is  proportional  to  x,  i.e. 

AH  « xA H,m.  (8) 


where  the  interchange  enthalpy.  A H,m.  is  a constant  characteristic  of  the 
considered  spinel  If  A S is  equal  to  the  configurational  entropy 
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X • Fraction  of  B loot  on  Tsfrahadral  Sites 
■ Degree  of  Inversion 

Figure  It.  The  dependence  of  the  confi- 
gurational entropy  in  spinels  on  the  degree 
of  inversion,  .v  (Reprinted  with  permission 
from  A Navrotsky  and  O.  J Kleppa. 
J.  Inorg.  Nucl.  ( hem..  29,  2701  (1967) 
Pergamon  Press  Ltd.) 


AG  = A H - TAS 

By  imposing  the  equilibrium  requirement  dAG/dx  = 0.  we  obtain 

_ . ,n 

RT  (1  - x)(2  - x)’ 


(10) 


(ID 


i.e.  a relation  between  x.  A//inl,  and  7 Navrotsky  and  Kleppa  (1967)  applied 
this  model  to  MgA)204.  and  made  a direct  determination  of  A/7  for  a natural 
spinel  by  means  of  oxide  melt  solution  calorimetry.  From  these  results  and 
other  information  they  concluded  that  for  MgAI204  at  1000  C,  A Hin,  * 10 
kcal/mole  and  x * 0. 1 . 
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More  recently  a modified  form  of  this  simple  disordering  model  was  adopted 
by  Holdaway  (1971)  in  order  to  resolve  some  of  the  internal  discrepancies 
between  the  experimental  observations  and  the  thermodynamic  data  for  the 
kyamte  sillimanite  andalusitc  (AI2Si05)  system.  Specifically.  Holdaway 
assumes  partial  cation  exchange  between  Al  and  Si  atoms  located  on  then  res- 
pective tetrahedral  sites  in  the  sillimanite  structure,  according  to  the  reaction 

AUi  + Sisj  = SiA1  + Als,  (12) 

His  principal  numerical  input  was  his  experimentally  measured  d P d 7 slope 
of  the  sillimanite  andalusite  phase  boundary,  which  gave  him  an  interchange 
enthalpy.  A of  about  14.75  kcal  mole.  With  this  interchange  enthalpy, 
one  would  predict  about  20  per  cent  disorder  in  sillimanite  at  1400  C. 

Quite  recently  oxide  melt  heat  of  solution  measurements  on  samples  of 
natural  sillimanite  heat-treated  at  1200  1700  C at  pressures  of  lb  23  kbar 
have  been  reported  by  Navrotsky.  Newton,  and  Kleppa  (1973).  The  results 
are  show  n in  Figure  12  and  indicate  a distinct  enthalpy  of  solution  decrement 
relative  to  unheated  sillimanite  of  about  1.3  kcal  mole  by  samples  run  at 
1 400  1550  C.  Analysis  of  these  results  in  terms  of  the  mentioned  disordering 
model  yields  a value  of  A Hml  =16+1  kcal,  mole,  in  good  agreement  with 
Holdaway's  value  derived  on  entirely  different  grounds.  At  temperatures  above 
1550  C.  larger  heat  of  solution  decrements  were  observed.  It  is  possible  that 
these  are  due  to  some  unknown  and  even  more  profound  disordering  process. 


600  800  1000  1200  1400  1600  1800 

Temperoture , °C 


figure  12  Heal  of  solution  data  on  quenched  samples  of  silltmanite  heat- 
treated  at  high  pressure  Solid  line  is  calculated  from  simple  Al  Si  tetrahedral 
disorder  Itq  |12||  with  lb  kcal  mole  tReprinted  with  permission  from 

A Navrotsky.  R ('  Newton  and  O J Kleppa.  (ieorhim.  Cosmochim.  Ana.  37. 

24V?  | IV73l  Pcrgamon  Press  Ltd.) 
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Id I M ineraloqieal  Reactions  and  Equilibria 

The  discussion  above  has  focused  on  only  one  of  the  serious  problems 
associated  with  the  application  of  thermochemical  methods  to  mineral  reactions 
and  equilibria,  i.e.  the  problem  of  structural  disorder.  Another  problem,  which 
we  have  not  discussed,  has  to  do  with  the  effect  of  chemical  impurities.  For 
example,  if  we  have  thermodynamic  data  for  pure,  ordered  phases  only,  how 
do  the  results  of  a thermodynamic  calculation  apply  to  a real  system  in  w hich 
significant  amounts  of  impurities  as  well  as  structural  disorder  may  be  present? 
Along  the  same  line  of  reasoning,  w hat  will  be  the  influence  on  the  considered 
equilibrium  of  more  extensive  solid  solution  formation  or  deviation  from 
stoichiometry,  which  at  times  may  vary  significantly  with  changes  in  tempera- 
ture and  pressure?  While  it  is  possible  to  give  partial  and  qualitative  answers 
to  some  of  these  questions  on  the  basis  of  chemical  knowledge  and  physical 
intuition,  more  quantitative  answers  generally  will  require  detailed  thermo- 
dynamic information  on  the  phases  which  actually  participate  in  the  considered 
reaction.  As  an  example  let  us  consider  the  hypothetical  reaction 

A + B = AB 

If  carried  out  in  the  laboratory,  this  reaction  may  well  be  associated  with 
solid  solubility  of  A in  B and  of  B in  A.  various  impurities  in  both  A.  B and  AB. 
and  perhaps  some  structural  disorder  and  solid  solubility  in  the  compound  AB 
Since  the  extent  of  these  effects  may  be  unknow  n or  difficult  to  assess,  it  may  be 
very  difficult  to  make  a realistic  quantitative  prediction  about  the  equilibrium. 

Let  us  assume  that  this  reaction  is  very  sluggish  and  can  be  studied  only  at 
very  high  temperatures  and  pressures.  Through  equations  (2)  and  (3)  this  will 
in  effect  furnish  a single  value  of  AG".  However,  the  equilibrium  work  does 
not  cover  a sufficient  range  in  temperature  to  establish  the  slope  of  the  P 7 
curve.  Therefore,  without  information  on  AS  or  on  AH,  we  are  unable  to  trace 
the  P I curve  to  lower  temperatures  and  pressures.  At  this  point  solution 
calorimetry,  if  carried  out  on  the  phases  which  actually  participate  in  the  reaction 
li  e.  on  the  reactants  A + B and  on  the  product  AB),  may  provide  information 
which  w ill  permit  us  to  calculate  AH.  and  hence  also  AS  and  the  P T slope 

The  procedure  which  has  been  outlined  here  for  the  hypothetical  reaction 
A + B = AB  is.  in  fact,  that  adopted  by  Newton  (1972)  in  his  work  on  the 
breakdown  of  cordierite  to  sapphinne  -t-  quart/  ( + enstatite '.’I.  More  recent 
and  more  extensive  work  by  Newton,  Charlu.  and  Kleppa  IW74)  has  led  to 
a much  more  detailed,  quantitative  understanding  of  this  reaction. 
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